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ABSTRACT 



Analytic expressions for distance-redshift relations which have been corrected 
for the effects of some inhomogeneities in the Friedmann-Robertson- Walker 
(FRW) mass density are used to illustrate the significance of inhomogeneities 
on a determination of go made by using Type la supernovae. The value of qo 
inferred from a given set of observations depends on the fractional amount of 
matter in inhomogeneities and is up to 50 percent larger than that obtained by 
using the standard Mattig m-z result for pure dust FRW models. 



Subject headings: cosmology: theory - large-scale structure of universe 
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1. Introduction 

Some years ago one of us (Kantowski 1969, henceforth referred to as Kl) used the 
"Swiss-cheese" models to study the effects of local inhomogeneities in the FRW mass density 
on the propagation of light through an otherwise homogeneous and isotropic universe. The 
Swiss-cheese models are standard pressure-free FRW models (sometimes called dust models) 
from which non-overlapping comoving spheres of zero-temperature material (called dust) 
have been removed and each replaced by a mass concentrated at its center. At a given epoch 
the model looks like Swiss-cheese full of perfectly round holes each with a seed (representing 
a spherical mass) at its center (see Fig. 1). If the holes are small compared to the radius of 
the FRW universe, the mass of the seed differs indistinguishably from the integrated mass 
of the dust missing from the sphere (see |K1| ). The comoving spheres expand/contract right 
along with the universe itself. The gravity field inside each expanding sphere is given by a 
Schwarzschild metric of concentrated mass. The importance of Swiss-cheese models is that 
they are explicit exact solutions to GR (with all boundary conditions being satisfied). In 
such models, effects of local inhomogeneities on observable relations can be calculated free 
of superposition approximations. 

The large scale structure of the universe, i.e., which pressure free FRW model we 
are actually living in, is clearly unaffected by how the clumps are distributed in number 
density or size. The FRW structure is fixed by giving the Hubble (Hq) and the deceleration 
(go) parameters. However, as we will show, observational determination of go is critically 
affected by the clumps. Most observations involve collecting radiation from a finite-sized 
distant source with a much smaller (~ point) receiver. The cross section of the radiation 
(the beam) ultimately received is consequently of the same order as the source's size. It is 
this size that is important in deciding what matter is to be ascribed to the seeds and what 
to the cheese. Matter whose size is small compared to the beam (e.g. an atom) is included 
in the cheese (the homogeneous dust) and matter whose size is large compared to the 
beam (e.g. a galaxy) is put into the seeds (the inhomogeneous clumps). The reader should 
remember that this is a model which is intended to mimic the optical effects of matter 
and not its actual physical distribution. For example atomic gas clouds are not physically 
confined away from galaxies as is the cheese; however, the optical effects of the gas should 
be independent of its location (see Section [2]]). In this model matter the same size as 
the beam clearly causes trouble and has to be dealt with separately for each observation. 
The above separation of mass based on size is not completely well defined because of the 
changing beam size from source to receiver. Star clusters for example would be considered 
small near a galactic sized source but large near us. 



In Section |2], optical propagation equations, originally derived in 



ikh 



, are restated. 



-4- 



These equations describe the propagation of a beam of electromagnetic radiation traversing 
a clumpy universe. The reader is directed to IK1] for details and to Sections 4.5.2 and 4.5.3 



of [Schneider et al. 1992| for a complete review of more recent work. In the Appendix it 
is concluded that the net effect of the seeds is negligible for z < 1 (except for uncommon 
lensing type situations) and in Section |3] it is shown that the resulting 6 = propagation 
equation for the cross sectional area of the beam is the associated Legendre equation. 
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The main result of this paper is the clear presentation of the corrected luminosity 
distance [expression ([18]) combined with (|i~4|)1 and the clear demonstration of the necessity 
of its use when evaluating q . In Section f|] the importance of taking into account the effects 
of inhomogeneities is demonstrated by showing how a recent determination of qo (made by 
applying the magnitude-redshift relation to SN1992bi) will be altered. 



2. Swiss-cheese Optics 



Observable relations such as magnitude-redshift or angular size-redshift are determined 
from expressions for the cross sectional area A of a light beam that propagates between 



source and observer. In [Kl| an integral-differential equation for the average area A as it 



propagates through a Swiss-cheese universe (see Fig. 1) was given as 

va' <a PD 

H — — — K>F > 



where 



y/A A 2 p F 



' - ! = -f< 1+ *>'< 1+ W*!' < 2 > 

(?) = iLH> !" A 2 {l+zfdv, (3) 
Pf Jo 



2 Before this manuscript could be completed and submitted and unbeknownst to us. 
equation ([□]) and its solutions ( |i~8| ) and (^T|) were published by [Seitz fc Schneider 1994 



They deserve credit for first recognizing (0) as the Associated Legendre equation and giving 
its general solution. They also show that the general solution reduces to previously known 
special solutions for appropriate parameter values. We give our version of the solution for 
completeness only. The conjectured extension of the validity of the 6 = equation beyond 
the Swiss-Cheese models is argued by |Schncider et al. 1992| in Sections 4.5.2 and 4.5.3. 
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-3g ^(l + ^) 5 , (4) 
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and where Weinberg sign conventions have been used ( [Misner et al. 19731) . 111 (HDj Pf is the 
usual FRW mass density (the local mass density of the cheese) and in (f|) Pf(0) is its current 
value. The other two mass densities, pr, and pi, are the components of p F made up of 
dust and inhomogeneities respectively, Pf = Pd + Pi- In the Swiss-cheese model pi is given 
by summing all clumped masses in a large volume (containing a representative number of 
clumps) and dividing by that volume. The density po is similarly computed by summing 
only the dust mass. The two driving terms for the area A in ([!]) are proportional to 1Z F 
and (£ 2 ) which represent average converging effects on the beam as it passes respectively 
through dust and near clumps. To arrive at ([!]) a large number of 'small' encounters with 
clumps had to be assumed. Any single encounter that causes a 'large' fractional change 
in the beam's observed cross sectional area was excluded. Such encounters cause lensing 
effects and must necessarily be accounted for by other means QBourassa fc Kantowski 1975| ). 
The v parameter is a particular affine parameter which has been expressed in (0) as a 
differential operator in redshift z. Equations (|2|), (^), and (|]) are slightly simpler than the 
corresponding equations in [Kl| because of the choice of the particular affine parameter made 
here (note that no loss of physical generality is made by such a choice). The electromagnetic 
wave starts from a point on the source at z = z s where v = 0, and expands at first as a 
small part of a spherical wave front, e.g. as a disk. As the spherical wave front encounters 
a clump the disk is distorted into an ellipse. The (£ 2 ) term in (TJ) is the effect on the wave 
front's area caused by repeated distortions of the clumps. The parameter b measures the 
strength of the distortions and has units of (length) -1 . Its value is defined by 

'(Gm/c 2 ) 1 / 3 '" 



H 2 ((Gm/c 2 )^) 
q ° c 2 ((Gm/c 2 )) 



(5) 



which makes b the order of the cube root of the Schwarzschild radius of a typical clump 
divided by its opaque radius squared (i.e., by a 2 ) and by the Hubble radius of the universe 
squared (i.e., by c 2 /Hi). 



The reader is directed to |K1| for details and assumptions made in all derivations, 
including the somewhat unphysical assumption that m/a 3 is the same for all clumps. If 
6 = 0, (|I|) is easily seen to be a linear homogeneous second-order equation in \J~A~ [ i.e., see 
(0)]; however, if b ^ 0, upon differentiation, (|I]) becomes a linear homogeneous third-order 
equation in A, 
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A'" - 4 



Pd 



Pf 



A' 



PD, 



Pi 



A = 0. 



(6) 



-K' F - —b a (l + z 
_Pf Pf 

The area A in (H) and @, when evaluated at the observer, is related to the luminosity 
distance D e by 5QD 2 = A\ (l + z s ) 2 or equivalently 59D e = y/A\ (l + z s ). The luminosity 
distance is defined in terms of the flux received T and the power 8V radiated into solid 
angle SQ by T = 5V/5QD 2 . The solid and plane angles are constants measured at 
the source (z = z s <-» v = v s = 0) for a wave front expanding into the future. The 
initial conditions for ([D) are VA\ S = and dy/~A/dv\ s = 86(1 + z s ) (or equivalently 
d\/A~/dz\ s = —58c/[H (l + z s ) 2 y/l + 2q z s \). The additional boundary condition needed for 
(0) comes from (|3|) and is d 2 \fA/dv 2 \ s = 0. These boundary conditions are a consequence 
of our choice of affine parameter v in (0), so chosen to simplify "apparent size" boundary 
conditions (see the next section). The general power series solution with b ^ is: 



VA\ /59 



Dt/(1 + z) = — 

-"0 




i + go 



z z + - 
2 



Pi i 
i + — go + % 

Pf 



Pi 
Pf 



-go 



c 3 pi 



—b 3 

30 H 3 p F 



z 4 + 0[z 5 



(7) 



where the redshift z is now the source's redshift. This result agrees with the Mattig result 
( IMattig 1958|) , 



D 



M 



rRoil + z) 



H,4 



(q -l) (y/l + 2q z - ij 



when pi = 0. 



If b can be neglected the difference between the standard Mattig result ([8p and (0) 
[or the exact result (|TH)] is due to the diminished converging effect of the mass remaining 
in the beam, i.e., the (pd/ Pf)7Zf term in (|I|). The standard result assumes all matter is 
uniformly distributed in and out of the beam. If some fraction pi/p F exists only out of the 
beam where it acts only through a negligible b 3 term, then a diminished fraction Pd/pf 
remains to converge the beam. In the Appendix we conclude that for a large majority of 
observations the b 3 term can be neglected. We now explicitly give the analytic solution of 
(D when 6 = 0. 



3. The Analytic Solution For A(z) when 6 = 
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The solution of ([I]) for 6 = and po = was given by Dyer fc Roeder 1972 and is the 
limiting case (pu — > 0) of result ( |T%D below. We first start by rewriting ([!]) when 6 = as 



(1 + z) 3 y /l + 2q z^(l + ^) 3 v /l + 2g ^ v /I(^) + ^3go(l + z) 5 /*^) 



0. 



(9) 



This equation^] can next be put into a recognizable self-adjoint form by changing the 
independent variable from z to £ and the dependent variable from J A(z) to P(C), 



(l-±- 

V 2g 0/ 

(i-CMO- 



i-C 2 



1 + 2 EE [ 

(l + zfyfAjz) EE (1-CM0- (10) 
For g < 1/2 we have £ 2 > 1 and for q > 1/2, is imaginary The resulting equation 



C 2 ) 



d 2 p 



^C 2 



-2C 



dp 



+ 



6 



Pi 



PF l - c 2 



p = o , 



(ii) 



has solutions which are arbitrary linear combinations of the associated Legendre functions 
of the first and second kind ( [Abramowitz fc Stegun 1964|) i.e., 



P(C)= Cl P 2 (C) + c 2 Q 2 (C) 



where 



l + 24 Pj /p F -l) /2. 



(12) 



(13) 



is the positive root of v{y + 1) = 6pj/p F . The range for v is from v — (usual pure dust 
FRW) where no clumps are present, to v = 2 where all mass is in the form of clumps and 
no dust is present. 

Boundary conditions can now be chosen appropriate for either luminosity distance Di as 
above or for apparent-size/angular distance D /_ . For apparent-size distances y/A = 58 D ^ , 
the boundary values are specified at the observer, P/|o = and dD^/dv\ = —1 (or 
equivalently dD/_/dz\o = c/Hq). For the 6 = case, the two types of solutions are easily 
related by the Wronskian of (§) (see |Etherington 1933] ), 



D e (z) = D L {z){l + zf. 



(14) 



3 See the Appendix A of Seitz & Schneider 1994 for the first presentation of this material 
including a comparison of all previously known limiting solutions. 
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Now switch to the more easily implemented apparent-size boundary conditions and hence 
require: 



(ciP*(C) + c 2 Ql(C) 

C\ 77 r C2 - 



d( 



d( 



z=0 



2=0 



o, 

c 

#o 



(15) 



These two equations can be used to find c\ and c 2 . The results then simplify by using the 
constancy of the Wronskian of (|TT|), 



to give 



d( dc (l-C 2 ) 



c 1 2 n2 / I 

c i = Tf — /i o 7 — ^V7 — I 7T7 w 7\Qu 



02 = ~JT vr^Wo Jy +W+ W)iy - l f 2y (tt^o) ■ (17) 

The resulting expression for apparent-size distance is 

J ^* ) = ivi- W + *K» + i)W(»-i) x <18) 



1 



1 + 2 



^2 f 1 \ p2 / Vl+2 g0 ^ _ 2 / 1 \ 2 / y/T+2^1 



This result is valid for all go > but some care has to be taken with the limit go ^ 1/2 as 
well as with the v — > 0, 1 limits. The z/ — > limit (pi = 0) is the Mattig result for FRW and 
is related to (^) by (|TJj). The v = 2 (i.e., the pr> = Dyer & Reoder 1972 result) is also 
obtained from (|I8D [see |5eitz fc Schneider 1994|| . 

Because the associated Legendre equation (|TT|) is a special case of the hypergeometric 
equation the above solution (|T8| ) can be written in terms of hypergeometric functions.^] This 
makes the analytic continuation to go > 1/2 transparent and gives a resulting expression 
for D i_ (or D$) whose constituent parts no longer have singularities. Computer routines 



x That differential equation ([!]) was of the hypergeometric type when 6 = was recognized 
by Pyer fc Roeder 1973| . 
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can then be used to evaluate these distances without having to avoid the singularities at 
q = 1/2, v = 0, and v = 1. The two identities needed are:[| 



p 2 m - r( ~^ ~ u ^ u+l F i v 1 % + 3 1 



2 



2^+1^1/2^2 _ i) z 1 2' 2' 2' C 

2T(± + z/)C +2 i/ _ 1 1 1 \ 

+ 7rV2r(z/ - 1)(C 2 - 1) 21 \2~ 1 '~2~2'2~ l/ ' , e)' ^ 



o*(c) = " 1/2 r(^ + 3)C- +1 p (u_l 3. J.' 

Vvig 2 -+ 1 r(z/+f)(C 2 -l) 2 2 ' 2 ' 2 'C 2 



7T 



1 / 2 r(z/ + 3)C^ 3 (C 2 -l) „ (v v 3 3 1 



"3T 2^1 k + 2, - + -; + -; - J , (20) 



2"+ 1 r(i/ + 1) V 2 2 2' 2' c 2 

giving the second form of the solution to equation @, 

c 1 

D L {q Q ,v-z) = — x (21) 

-no + 2) 

l + 2q z) 1+ ^ 2 „ (v , n 1/ 3 3 n n \ „ { v 1 f 11 1 - 2g \ 



2 J 



Fi - + 2, - + -; v + -; 1 - 2g Q 2 Fi -- - 1, -- - -; - - v 



1 + z) 2 \2 '2 2' 2' "V V 2 2 2' 2 ' 1 + 2q z i 

1 / v is 1 1 \ (v v 3 3 l-2o 

2 Fi(---l,----;--i/;l-2g UFi -+2,- + -;i/+ ' 



(1 + 2g z) 3 / 2 +<V 2 z 'I 2 '2 2' 2 ' w / V 2 2 2' 2 , l + 2q z 

This form is well behaved as a function of z/ as well as in the neighborhood of 90 = 1/2, 
and is needed beyond the circle of convergence (|£| = 1) of the hypergeometric series 
2-^1(0, b; c; £) only on the negative real axis £ = x < — 1 when q > 1 (where it is easily 
continued). In Fig. 2, D /_ is plotted as a function of z for three values of v and three 
values of q and in Fig. 3 the related magnitude-redshift curves are plotted. The reader 
should recall that apparent (bolometric) magnitudes m are related to absolute (bolometric) 
magnitudes M and luminosity distances Di [see (0)] by 

m = M + 5 \og 10 (D e /10pc) = 5 log 10 (-^A>) + constant, (22) 

which when evaluated as a function of redshift is called the m-z (magnitude-redshift) 
relation. As pointed out in the next section the constant = M — 5 log 10 (-^-10pc) turns 
out to be observationally independent of the Hubble parameter. In the next section we 
investigate the sensitivity of this relation to the dumpiness parameter v. 



Note that in 8.1.5 of |Abramowitz fc Stegun 1964| a it 1//2 is missing 
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4. Application 



We conclude by illustrating the effect of inhomogeneities on the determination 
of qo when Type la supernovae (SNe la) are used in the Hubble diagram. Such a 
determination is expected to be feasible in the near future. Attempts to use galaxies 
or quasars in this way are foiled by the intrinsic dispersion in their luminosities and 
by uncertainties associated with their evolution. However, normal SNe la (those that 
have nonpeculiar spectra, light-curve shapes, and broad band colors) are luminous point 
sources that form a well defined ridge line in a Hubble diagram (for recent reviews of 
SNe la see [Branch fc Khokhlov 1994| and references therein). Even just the application 



of a simple B-V color criterion defines a subsample of SNe la that has mean absolute 
magnitudes Mb — My — —19.8 ± 0.1 + 51og(Ho/50) and dispersions about the means 
cr(M B ) ^ cr(M v ) ^ 0.2 magnitudes QVaughan et al. 19951 , [Hamuy et al. 1995) , |Vaughan~fcl 



[Branch 1995| ). With sufficient spectroscopic or photometric information it may be possible 
to control the individual SN la relative magnitudes to 0.1 or even better ( Fisher et al. 



1995| , [Hamuy et al. 199"o] , [Riess, Press, fc Kirshner 1995| , [l'ammann fc Sandage 1995] ). The 



value of H , and thus the calibration of the mean SN la absolute magnitude, does not enter 
the determination of go from the Hubble diagram. 

In the standard SN la model, a mass accreting white dwarf approaches the 
Chandraskhar limit, ignites degenerate carbon near its center, undergoes a thermonuclear 
instability, and disrupts completely. The light curve is powered by the radioactive decay 
of freshly synthesized 56 Ni through 56 Co to 56 Fe. Thus, to first order, the peak absolute 
magnitude of SNe la is not expected to vary with cosmic time. Observationally, some 
evidence for a dependence of the properties of SNe la on the age of the stellar population 
at their sites has been found ([Branch fc van den Bergh 1993| , [Hamuy et al. 1995f) . This 
constitutes evidence for mild evolution of the SN la population with cosmic time, but it 
implies that SNe la at high redshift, when the universe was younger, ought to have an even 
smaller dispersion among their properties than do the nearer SNe la that constitute the 
bulk of the present observational sample. In any case, with sufficient spectroscopic and 
photometric information it should be possible to identify counterparts of the high redshift 
SNe la among the nearer SNe la. 

Recently it has been demonstrated that SNe la at cosmologically significant redshifts 
can be discovered and measured. |Perlmutter et al. 1994| has announced the discovery of SN 



1994F, a SN la at z = 0.354; SN 1994G, a probable SN la at z = 0.425; and SN 1994H, a 



probable SN la at z = 0.373. Perlmutter et al. 1995| present a light curve of SN 1992bi, a 



probable SN la at z=0.458, which they use with the standard Mattig relation to make a 
preliminary estimate of the value of go- Recently additional remote SNe la and probable 
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SNe la have been discovered (C. Pennypacker, personal communication). 

Fig. 4, which is an enlarged portion of Fig. 3 near the redshift z = 0.458 of SN 
1992bi, shows that clumping can have a significant effect on the determination of go a t this 
redshift. In Fig. 5, the quantity 51ogio (H /c)Di(q ,i / ; z = 0.458) [see (p2|) 1 is plotted as a 
function of go for three values of the clumping parameter v. For example, if a given set of 
observations gives go = 0.1 (|Perlmutter et al. 1995] ) when the Mattig relation is used, the 



v = 2 (completely clumpy, pi = pp) relation would give g = 0.15; that is, if the universe is 
completely clumpy then neglecting to take dumpiness into account causes an underestimate 
of g by about 50 percent. The difference between g = 0.1 and g = 0.15 is significant 
quantitatively, but not qualitatively. On the other hand, if the observations were such that 
the Mattig relation gave go = 0.34, then the v = 2 relation would give go = 0.51. This 
difference is more than merely quantitative — it is the difference between an open and a 
closed universe. 

We find that it is just when redshifts are large enough for the determination of go 
that the effects of clumps become significant. Figs. 4 and 5 show that if g were known, 
distinguishing between a completely clumpy universe and a universe with no clumps would 
require that the average apparent magnitude of a set of remote SNe la be compared to that 
of nearby SNe la to an accuracy 5m ~ 0.05 magnitude. Such statistical accuracy should 
soon be available at the rate SNe la are being discovered. In practice, though, the m-z 
relation only constrains the go and v parameters in combination. Independent information 
must be used to determine their individual values. 



This work was supported by the Department of Energy, the Southern Association 
for High Energy Physics (SAHEP) funded by the Texas National Research Laboratory 
Commission (TNRLC), and by NSF grant AST 91-15061. 



A. Can the parameter b be neglected? 



Based on the series solution (0) , [Kl] concluded that the b 3 term appearing in (HD and (|) 
would be insignificant for observations of large objects such as galaxies. Further arguments 
are now given to reassert that the probability of the b 3 term affecting an observation is 
rather small. A priori, observing a significant effect of b on the magnitude-redshift relation 
should be about as likely as a single observation resulting in a gravitational lens effect. 
When either lensing or a significant distortion of a light beam occurs, the cause is the 



-12- 



same conformal curvature term in the gravity field. Rauch (1991) has come to a similar 
conclusion about the unlikely probability of observing gravitationally lensed SNe la at 
small z as we do here. In order to estimate the significance of focusing by external mass 
clumps on a passing beam of light, i.e., to estimate the relative significance of b 3 in (|B|), 
a comparison of the magnitudes of the two terms that appear as coefficients of A can be 
made. Even if such a comparison implies that the b 3 term is significant compared to the 1Z' 
term, they both might be unimportant compared to the z development caused by the A'" 
or A' terms. The ratio of the two A terms is 

Pi b 3 {l + zf = pj_ ( b 3 \ 1 fA1 , 

Pd W p D \l5q (H /c) 3 ) (1 + z)y/l + 2q z ' 1 ' 

The unitless combination 



b 3 _ 4 c ((Gm/c 2 ) 4 / 3 ) 
0= 15go(#o/c) 3 ~ 5H ((Gm/c 2 )) 
is clearly the parameter determined by the structure of the inhomogeneities, which affects 
the relative importance. In the series solution (|?p, Bq first appears in the z A term, already 
indicating that it is probably not important at small redshifts. 

Dyer & Roeder 1974 numerically integrated @ for numerous values of Pi/pd and 
B = 15Bq. From their graphs we can conclude that as long as Bq ~ 1 and observations are 
at redshifts z < 1 then the b 3 term is not significant. In Fig. 6 the recomputed numerical 
solutions of (|) are given for values of B = 0, 1, and 5, go = 0.1 and 0.5, and Pi/pf = 
and 1, confirming the above conclusion. For example for g = 0.5 the magnitude of a SN 
la at z ~ 0.5 in the standard dust filled FRW universe would appear 0.075 magnitudes 
brighter than in the Bq = clumpy universe. If, in the clumpy universe, Bq = 1 or Bq = 5 
the magnitude would be brighter by 0.01 and 0.05 magnitudes respectively. For g = 0.1 
all changes are about 1/4 of those for go — 0.5 . The conclusion is that for Bq ~ 1 (and 
smaller) the b 3 term makes only a small increase in intensity (~ 1/7 or less) compared to 
the decrease in intensity due to the missing matter in the light beam [y = 2 rather than 
v — 0). As to what Bq values are known to exist: for globular clusters, Bq ~ 10~ 2 ; for 
galaxies and galaxy clusters, Bq ~ 1; for stars, Bq ~ 10 9 ; and for black holes at the centers 
of spiral galaxies, Bq ~ 10 25 . 

What is now shown is that the probability of seeing the effects of a Bq larger than 
~ 1 is quite small (for these limited redshifts). If it is assumed that there is a randomly 
distributed population of essentially point inhomogeneities of number density n = uqRI/R 3 
between us and an observed object at redshift z, the probability that the observed light 
beam comes no closer than a distance a' to the center of any of the intervening masses can 



(Gm/c ) ' 



(A2) 
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be computed by evaluating 



P(r) 



cxp 



rma , Rdr 
VI - kr 2 



(A3) 



Assuming the clumps are points will maximize their possible effects. The parameter r used 
here is the FRW comoving radial coordinate centered on the source and R is the FRW scale 
factor. This result follows from the simple assumption that —AP = P(r) — P(r + Ar), the 
probability a beam make it to r having had no impact parameter less than a' and then 
having an impact parameter less than a' within the next Ar, is given by 



AP = P(r) 



a 



x n Airr 2 ArR 3 



(A4) 



47rr 2 i? 2 \/l - kr 2 

This is a simple ratio of the sum of impact disk areas in a spherical shell of thickness 
RAr/y/1 — kr 2 to the total area on the sky, at comoving distance r. After changing from 



co-moving coordinate r to redshift z, 



can be easily integrated to give 



P(z) 



exp 



7ra' 2 noC 



f(z,qo) 



where 



f(z,Qo) 




1 + 2q z) 3 / 2 - 1 



l-2q ) (l + 2q z 



,1/2 



(A5) 



(A6) 



If it is assumed that all members of this population of inhomogeneities have the same mass 



m, 



5|) can be rewritten as a function of the unitless parameter B' of ( |A2| ) as 



P(z) = exp 



3<?o pi 
5B' qPf 



f(z,qo) 



(A7) 



P(z) is interpreted as the probability that all impact parameters (a values as the light beam 
interacts with inhomogeneities on its path from a source at z to us) have values such that 
AGm I \5H ca 2 ) = B is less than B' . In Fig. 7, P(z) is plotted for several values of q and 
(pi I ' Pf)B'q. For a given B' Q the maximum probability occurs when pr/ P F — 1, i-e., when all 
matter is in point clumps. From Fig. 7 it is clear that for redshifts z < 1, the probability of 
having a value of B < 1 is high. For such values of B it has already been concluded that 
the b term in (|^) has a rather small effect. What is expected to happen is that an occasional 
observation of a SN la will occur where the intensity is increased because it happens to be 
closely aligned with a dark and unseen interveining galaxy. According to Fig. 7 this can't 
happen frequently. However, when it does, the SN will either sit far from the m-z curve 
where it can be excluded, or sit close to the m-z curve and simply increase the error in the 
estimated value of go- 
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Fig. 1. — Radiation beam of cross section A propagating through a Swiss-cheese universe 
from distant source to observer. 



Fig. 2. — Apparent-size distance -^-D<(go? v\ z) as a function of redshift z for three values 
of go and three values of v. 



Fig. 3. — Magnitude-redshift relation, 5 log 10 ^-D^qo, v \ z), as a function of redshift z for 
three values of q and three values of v. 



Fig. 4. — Magnitude-redshift relation, 5 log 10 ^Dt(q , v\ z) , as a function of redshift z for 
three values of q and three values of v. 



Fig. 5. — Magnitude-redshift relation, 5 log 10 — De(qo, u; z = 0.458), as a function of the 
deceleration parameter q for three values of v. 



Fig. 6. — Magnitude-redshift relation, 5 log 10 — Dt(q , u, B ; z), as a function of redshift z 
near 0.458 for two values of q and u, and three values of B . 



Fig. 7. — Probability of B < B' in a universe of given Pi/pd 




Fig. 1. — Radiation beam of cross section A propagating through 
a swisscheese universe from distant source to observer. 
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Fig. 2. — Apparent-size distance -^-D/_(qQ, u; z) as a function of 
redshift z for three values of go an d three values of v. 
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Fig. 3. — Magnitude-redshift relation, 51og 10 -^-D^qo, v\ 2), 
as a function of redshift z for three values of go an d three values of v. 

3 




Fig. 4. — Magnitude-redshift relation, 51og 10 -^-D^qo, v\ 2), 
as a function of redshift z for three values of go an d three values of v. 
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Fig. 5. — Magnitude-redshift relation, 51og 10 -^-D^qo, v\ z = 0.458), 
as a function of the deceleration parameter go f° r three values of v. 
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Fig. 6. — Magnitude-redshift relation, 51og 10 -^-D^qo, v, Bq; 2), as a function 
of redshift z near 0.458 for two values of go an d v, and three values of Bq. 



6 




0~ 0.5 1~~ 1^5 ~2 



Redshift z 



Fig. 7. — Probability of Bq < B' in a universe of given pi/po 
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